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An exact formalism for the relativistic version of Landau theory of Fermi liquid in presence of 
strong quantizing magnetic field is developed. Both direct and exchange type interactions with 
scalar and vector coupling cases are considered. 



1. INTRODUCTION 



o ■ 
o : 

The study of dense stellar matter under ultra-strong magnetic field has gotten a new life after the discovery of a 
5^ ' few magnetars [1-4]. These stellar objects are believed to be strongly magnetized young neutron stars. The surface 
magnetic fields arc observed to be > 10 15 G. Then it is quite possible that the field at the core region may go up to 
10 18 G. The exact source of this strong magnetic field is not known. These objects are also supposed to be the possible 
sources of anomalous X-ray and soft gamma emissions (AXP and SGR). If the magnetic field is really so strong, in 
particular at the core region, they must affect most of the important physical properties of such stellar objects. The 
elementary processes, e.g., weak and electromagnetic processes taking place at the core region should also change 
. significantly. 

H ' The strong magnetic field affects the equation of state of dense neutron star matter and as a consequence the 
gross-properties of neutron stars [5-8], e.g., mass-radius relation, moment of inertia, rotational frequency etc. should 
change significantly. In the case of compact neutron stars, the phase transition from neutron matter to quark matter 
at the core region is also affected by strong quantizing magnetic field. It has been shown that a first order phase 
transition initiated by the nucleation of quark matter droplets is absolutely forbidden if the magnetic field ~ 10 15 G 
at the core region [9,10]. However a second order phase transition is allowed, provided the magnetic field strength 
< 10 20 G. This is of course too high to achieve at the core region. 

The elementary processes, in particular, the weak and the electromagnetic processes taking place at the core region 
of a neutron star are strongly affected by such ultra-strong magnetic field [11,12]. Since the cooling of neutron stars 
are mainly controlled by neutrino/anti-neutrino emission, the presence of strong quantizing magnetic field should 
^"i , affect the thermal history of strongly magnetized neutron stars. Further, the electrical conductivity of neutron star 
■ matter which directly controls the evolution of neutron star magnetic field will also change significantly. 

In another kind of work, the stability of such strongly magnetized rotating objects are studied. It has been observed 
^ . from the detailed general relativistic calculation that there are possibility of some form of geometrical deformation 
V^J ' of these objects from their usual spherical shapes [13-15]. In the extreme case such objects may either become black 
strings or black disks. In the non-extreme case, however, it is also possible to detect gravity waves from these rotating 
C3 ' objects. 

In a recent study of microscopic model of dense neutron star matter, we have observed that if most of the electrons 
occupy the zeroth Landau level, with spin anti-parallel to the direction of magnetic field and very few of them are with 
spin along the direction of magnetic field and Landau quantum number > 0, then either such strongly magnetized 
system can not exist or such a strong magnetic field can not be realized at the core region of a neutron star [16]. We 
have shown in that work that since the electrical conductivity of the medium becomes extremely low in presence of 
ultra-strong magnetic field, the magnetic field at the core region must therefore decay very rapidly. Hence we have 
argued that the second possibility is more feasible, i.e., strong magnetic field can not exist at the core of a magnetar. 

So far most of the calculations on the equation of states of dense stellar matter in presence of strong magnetic 
field are either based on some kind of mean field approximation or non-relativistic potential model [6,17]. In this 
paper we shall derive an exact formalism of relativistic version of Landau theory of Fermi liquid in presence of strong 
quantizing magnetic field and obtain the quasi-particle energy for both scalar and vector coupling cases. We shall 
consider a typical many body fermionic system in presence of strong magnetic field interacting via some kind of scalar 
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and vector bosons exchange. Therefore from this investigation it is very easy to obtain equation of state of a real 
fcrmion system in presence of strong magnetic field, e.g., high density electron gas (with photon exchange), dense 
neutron star matter (with a — lo — p meson exchange) or even dense quark matter produced at the core region as 
a result of quark- hadron phase transition (exchange of color gluons). Not only that, one can study the magnetic as 
well as the non-equilibrium properties of such dense fermionic system as mentioned above. The relativistic version of 
this model without magnetic field was developed long ago by Baym and Chin [18]. So far most of the scattering or 
reaction or decay rates/cross-sections are evaluated cither with very low magnetic field approximation or in presence 
of ultra-strong magnetic field. To the best of our knowledge the present article is the first attempt to develop an exact 
theoretical formalism to obtain the two-body scattering matrix element (both direct and exchange) for all possible 
strengths of external magnetic field. In a completely different type of work we have used this formalism to obtain 
the rates of weak and electro-magnetic processes. We have investigated the effect of strong magnetic field on the 
emissivities and mean free paths of neutrinos / anti-neutrinos and electrical conductivity of dense electron gas in 
neutron star matter. These quantities play significant role in thermal evolution and the evolution of magnetic fields 
of neutron stars respectively. 

Now similar to the conventional non-relativistic case, the relativistic version also deals with normal Fermi liquid 
and applicable only for the low-lying excited states of the system, which are made of superpositions of quasi-particle 
excitations close to the Fermi surface. The relativistic version is therefore applicable to high density electron gas, 
dense neutron star matter and also to dense quark matter which may present at the core region of compact neutron 
stars. For the applicability of this model the temperature of the systems are therefore assumed to be low enough 
compared to the chemical potential of the constituents. 

The paper is organized in the following manner. In the next section we shall develop an exact formalism for the 
relativistic version of Landau theory of Fermi liquid considering a typical many body fermionic system interacting 
via some kind of scalar and vector boson exchange field. In the last section we have discussed the results and future 
prospects of the work. 



We have considered a dense relativistic fermionic system at zero temperature (T << [if, the chemical potential of 
the fermions) and interacting via scalar and vector boson fields represented by <f> and V M respectively with masses £ 
and r\ respectively. The corresponding coupling constants with the fermionic field are g s and g v respectively. Then 
the basic Lagrangian density is given by 



Where = <9 M + iqA^ with q is the magnitude of electric charge, m is the rest mass of the fermions and = 
(0, 0, xB m , 0), the electro magnetic field vector. Here the gauge choice is such that the constant external magnetic field 
B m is along Z-axis, C v , C s and £ cm are the free Lagrangian densities for the vector, scalar and the electromagnetic 
fields respectively. Therefore, in this particular case, £ om is the Lagrangian density corresponding to the constant 
magnetic field B m . We now adiabatically turn on the interaction between the fermions. Therefore, the normal Fermi 
liquid system will evolve from the non-interacting ground state into the interacting ground state and there will be a 
one-to-one correspondence between the bare particle states of the original system and the dressed or quasi-particle 
states of the interacting system. 

Next to obtain / the Landau Fermi liquid interaction, which is also called the quasi-particle interaction function, 
we have considered both direct and exchange type interaction. The exchanged bosons are either scalar or vector 
type. The Landau Fermi liquid interaction / is related to the two-particle forward scattering amplitude. Therefore, 
it is essential to compute transition matrix element for two-particle direct as well as exchange interaction forward 
scattering via the exchanged scalar and vector bosons. We further assume that the constant magnetic field acts as a 
background field and is present throughout the system. We further assume that the strength of magnetic field is such 
that in the relativistic case the Landau levels are populated for the charged fermions, i.e., B m > Bm^\ the quantum 
mechanical limit, given by qB$^ — to 2 . 

Then the modified form of spinor solutions of charged fermions in presence of strong quantizing magnetic field is 
given by 
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where the symbols t and I indicates up and down spin states. 
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with H v the well known Hermite polynomial of order u, E v = (p 2 z + m 2 + 2vqB m ) 1 / 2 , the single particle energy 
eigen value with v = 0, 1, 2, the Landau quantum numbers and L y , L z are length scales along Y and Z directions 
respectively. 

In presence of strong quantizing magnetic field along Z-axis (as specified by the choice of gauge) the momentum in 
the orthogonal plane gets quantized and is given by p±^ = {2vqB m ) 1 / 2 , whereas, the component along Z-axis varies 
continuously from — oo to oo. Further the phase space volume integral in momentum space in this condition is given 
by 



(2 



W I d3pf(p) = (^rJ dPzd2p±m M /_7 dpJi ^ p 
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(6) 



here we have assumed c = h — and pf is the Fermi momentum. 

To compute the Landau Fermi liquid interaction function / from two particle forward scattering matrix, we first 
derive from the first principle the general form of two fermion transition matrix element. We calculate for both scalar 
and vector coupling with direct and exchange type interactions. Therefore, we have the two particle transition matrix 
corresponding to scalar and vector boson exchange interactions 

(7) 
(8) 



T fi = i J j° fi (x)ct>(x)d 4 x 
T fi = i J j^(x)V»(x)d 4 x 



where the fermion four-current is given by jf^x) = gk'ipf{x)j ,1 tpi(x) = (p(x),j(x)), with the spinors given by eqns.(3) 
and (4) and k = s and k = v corresponding to scalar and vector coupling cases respectively. The scalar and vector 
fields can very easily be obtained from the well known Greens function solutions of Klein-Gordon equations 



and are given by 



<f>(x) 
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(d 2 + m 2 )cb,V^ p,f 
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(9) 
(10) 



where Q is the transferred four momentum (Q^ = (qo,q x ,q y ,q z )). Let us first compute the transition matrix element 
for scalar interaction case. Substituting for the Fermi spinors (eqns.(3) and (4)), we have 
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p(x') = -r f r exp[-i{(E„ 1 -E V2 )- ( Ply - k ly ) - ( Plz - k lz )}][p] (11) 

where 

W = \^u{x\v2M)l u(x',v u p{) = ^{h{x' ,vi,u 2 ,PiM)l°) (12) 

spin 

Then we have after putting for p(x') and integrating over t', y' and z', the direct part of two-particle transition matrix 
element with scalar coupling interaction 



l fl - 1 



CXp[-i{(£„ ; - E„> 2 ) - (p 2y - k 2y ) - (P2z - k 2z )}} (13) 

Integrating over t, y, z and q°, q y , q z , we get 

T)f s) = - l (2nfS(E)S(p y )S(p z ) J ^ J ^[ p ( x ) Mx ')] cxp[ - lqx ( x - x ')]^J-_dxdx' (14) 

where the (5-functions indicate the abbreviated form of energy, y and z components of momentum conservation. In the 
derivation of eqn.(14) and the computation of subsequent matrix elements for direct case with vector boson exchange, 
we have assumed the two body process in the form 

l(pi)+2(p 2 )-fl / (fci) + 2 / (fe2) (15) 

In this case the forward scattering amplitude can be obtained if we put E Vl = E v , 2 = E v (say), E v ^ = E v ^ = E v i 
(say), pi y = ki y = p y (say), p 2y = k 2y = p' y (say), p\ z = k\ z = p z (say) and p 2z = k 2z = p' z (say). Then after 
evaluating the contour integral over q x , given by 

T _ [ +oc dg x 1 _ exp[-C I x - |] 

^-J^ 2nQ*-e- 2C (ibj 

and substituting in two-body forward scattering matrix element given by eqn.(14), we have the Landau Fermi liquid 
interaction function corresponding to direct term for scalar coupling case 

= / ex PK \x~x' \}[p(x)}[p(x>)}dxdx> (17) 
After evaluating the quantities within [ ] (see Appendix A) wc finally get 

•>■) r r r , \ / ttt ~T — r / cxp[-C \x-x' \]dxdx' 
2,2C,E v E v i{E v + m)(E u > + m) J 

{2E V {E V + m){ll. py (x) + ll_ x . Vy (x)) + 2vqB m (ll Vy (x) - I 2 _ 1;Py (x))} 

{2E V ,(E V , + m){ll,. ply {x') + #_ 1;p , (a/)) + 2v'qB m {ll,. p , y {x') - #_ 1;p , (*'))} (18) 
Then the quasi-particle energy is given by 



As) _ gl_ 

JDii 



£ v (p z ) = E„(p z ) + ——2 ^2 (2-5„> ) / / / dxdx'dp y >dpz>fvlXP,p'\v,v') 

V^ 71 / j,' = o Jx,x'= — oo Jpyi = — oo Jp_,,=~ pf 

= E v {p z ) + ^E% T XPz) (19) 



where [f max ] = (p"j — m 2 )/(2qB m ), the integer part of right hand side. In the case of ultra-strong magnetic field, 
v = v' = 0, then 
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q 2 

eo(p z ) = E Q ( Pz ) + |fn/ (20) 

where 

n f = ^TPf ( 21 ) 

the number density of fermions. The form of the result given in eqn.(20) is identical with that of zero field case [18]. 
The energy density is then given by 



qB 

e = 



|wj fPf 

(2~M / d Pz e v ( Pz ) (22) 



4tt 2 

i/=0 J -Pf 

Let us now compute the exchange part for scalar coupling case. In this case we substitute E Vl = _E, y < = E v , 
E^ = E U2 = E„>, piy = k 2y = p y , piz = k 2z = p z , Ply = ki y = p' y , p2 Z = ki z = p' z , to obtain the two-body forward 
scattering amplitude. In this particular example, the basic process is 

l(pi) + 2( P2 )-fl / (fe2) + 2'(fei) (23) 

and the q x contour integral is given by 

cxp[- K | x-x'\] 

where K « + q\ + C 2 ) 1 ^ 2 and % = (pt — p'A, with i = y and z. Then evaluating [p{x)p{x')] (see Appendix A), we 
get 



.,2 



= / dxdx 



,exp[— K | x — x' |] 

— — — i n.^nX 

2,2E V E V , J K 
\{E V E V , + PzP ' z + m 2 )Ti(AA') - p±.j? ± Tr(BB')] (25) 

where 

Tr(AA') = 2{I v {x)I v {x')I v ,{x)I vl {x l ) + I v - 1 {x)I v -i(x')I v .-i(x)I v .- 1 (x')) (26) 

and 

Tt(BB') = 2{I v _ 1 {x)I v {x')I vl . 1 {x)I v ,{x') + J„ (a/) (x) J,x_i (a;')) (27) 

Then as before 
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For i/ = i/ = 0, 



A25« - *L T ^ exp ^ (^+^+^ 2 ) 
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Pi 
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erfc( " 1/2 seed ) so. -fr/fl (29) 



where u' w| p z — p' z |. In the zero field case 



f (s) = 9 2 s m 2 ~p.p' 

Jcx 4E Q (p)E Q (p>)(p'-p) 2 + e [ ' 

We next consider the vector coupling case. The current four vector in the direct case is given by 
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n*) - 9v e -^tmi - ^) zpi - ^) - - ^ P1> fcl)] ( 3i) 

J-'yJ-'z 

where 

= ^E^" u= ^(A(z,^2,Pi, W) (32) 

spin 

Then it is very easy to show that the Landau Fermi liquid interaction function is given by 

fmr = QO p P / exp[-r? | x - x' \][{E V E V , - PzPz ,)TrA TrA' + ± .0 ± TrB TrB'] (33) 

where we have used the same kind of substitutions as has been done for the scalar case and further the relation, 

Tr(A(x, v, p)y ) = A- {TrA + TrB)p» (34) 

(see Appendix B) is used to obtain this expression. 

The exchange term corresponding to vector coupling case is similarly given by (see Appendix B for brief derivation) 

^ = 32ifk; / dxdx ' CM ~ Kl K X ~ X ' 'W iV - E » E »> + 3™ 2 )Tr(AA') - p ± .j? ± Tr{BB')] (35) 
Then energy density of the system is given by (including all the terms) 
qBr, 



4tt 2 



W [ p+pf r+p f | 

VJ (2 - 5,0) { / dp^fo) + ^ / dp^E^ipz) \ (36) 



where k = s or u and / = Dir or ex. The kinetic pressure is then given by 

P = m - e (37) 
The chemical potential or the Fermi energy is given by 

H = e v (p z = Pf ) = E v {p f ) + AE^pf) 

k,l 

= Mo + (38) 

Therefore knowing the modified form of single particle energy or the quasi-particle energy, it is possible to obtain all 
the thermodynamic quantities of the system from standard definitions. 

To obtain equation of state of dense neutron star matter or dense electron gas, since both direct and exchange terms 
contribute in the case ofe + e^e + e scattering, it is necessary to screen e — e and e — p interactions in the direct 
case (in the case of e — p scattering, of course, only the direct term exists) . To get the screening mass of electron, 
we proceed in the following manner. We assume protons as the positively charged background. Then assuming local 
charge fluctuation we have 

n + = rip + n+ 

n e = Uq + n' e (39) 

In the case positive lattice background, we have to replace n by Zn n , where Z is the atomic number of the lattice 
ions. Then in the case of proton background, the well known Poisson equation is given by 

V 2 ?A = -4:irnfe + 4irn' e e (40) 

where ip is the electrostatic potential produced in the system because of local charge fluctuation. Then we have from 
eqn.(21) after substituting for pf in terms of chemical potential [if and electrostatic potential ip, 

+ . eB m ^ fi f eip 



G 



The last term is the perturbed part (n' e ) as mentioned in the second line of eqn.(39). Then we have 

„2r „. i 



V 2 V - -47rn+ e + 26 Bmflf £(2 - 5 v0 )- 



(jij -m 2 -2veB m y/ 2 
which may be re-written in the form 

[V 2 - fc s 2 c ] V = -Annfe 
where k sc is the screening mass and is given by 

1/2 



7T ' * 
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(^ - m 2 - 2i/eB m )V2 



(42) 



(43) 



(44) 



and the screening length rr> = l/fc S c 

In the case of quark matter, however, only the exchange term contributes, we therefore do not need quark-quark 
screening. Since colorless gluons do not exist, the direct term has no significance. 

3. CONCLUSION 

In this article we have developed an exact formalism to obtain two-body scattering matrix element for an wide range 
of magnetic field strength. Although the result has been used to formulate an exact relativistic version of Landau 
theory of Fermi liquid, is also applicable to study weak and electromagnetic processes in presence of strong magnetic 
field. 

The formalism developed in this paper to compute Landau Fermi liquid interaction / can also be used to study 
magnetic properties of fermionic system. The calculations (both theoretical and numerical) are of course quite 
complicated in the relativistic region (we shall report the result in some future publication). The other interesting 
property- the transport theory of normal Fermi liquid can very easily be studied with the help of present formalism 
and obtain various kinetic coefficients along with their dependences on the strength of magnetic field. 



APPENDIX A 



In the direct case we need 



A(u, x, x', k y ) — u(u, x, k)u(u, x', k) 



spin 



Substituting for the up and down spin solutions of Dirac equation (eqns.(3) and (4)), we get 

A = ^r{Ak^{ji = and z) + mA + Bk^(/j, = y and p y = p±)) 
The matrices A and B are given by 
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where the primes indicate the functions of x' . 
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Eqn.(A2) is an entirely new result and to the best of our knowledge it has not been reported earlier. Further, this 
result plays the key role in all kind of calculations related with electromagnetic and weak interactions in presence of 
strong quantizing magnetic field. 

Since 7 matrices are traceless and both A and B matrices are diagonal with identical blocks, it is very easy to 
evaluate the traces of the product of 7- matrices multiplied with any number of A and/or B, from any side with any 
order e.g., 

Tr( 7 ^ 7 M lJ 4 2 ..B 1 B 2 ..) = Tr(A 1 A 2 ..B 1 B 2 ..)g^, (A5) 

Tr( r f^ 1 °A 1 A 2 ..B 1 B 2 ..) = Tr^A^.B^..)^ g^ - g» x g v ° + g^g""), (A6) 

Tr(product of odd 7s with A and/or B) = etc. The other interesting aspects of A and B matrices are: 

i) k 1 ^Tr{A 1 A 2 ) = {E X E 2 - k lz k 2z )Tr(A 1 A 2 ) 

ii) k^TriB^) = fcix-fexlir^iSa) 

hi) k^k^Tr^Bi) = fc lAl fc^Tr(£M 2 ) = 

iv) pi tl k 1 ^p 2v k 2v Tv{A 1 B 2 ) ^ = (E Vl E u ^ —pi z ki z )p 2 ±.k 2 ±Tr(AiB 2 ) These set of relations are also totally new results 
and have not been reported before. 

To compute the direct part for scalar coupling case, we need 

Tr(u(a;, v, p)u(x, v, p)j a ) = Tr(A(.x, v, p)j°) (A7) 

Using the expression for A, it very easy to show that the above trace is given by 

2EA E v + m) VMEu + m)(ll py (x) + ll_ V;py (x)) + 2vqB{ll Py (x) I*_ 1;Py (x))} (A8) 

To compute the exchange part for this case, we need 

Tv[(u(x, i/,p)7°«(i/, x,p'))(u(x', is',p'h°u(v, x',p))} = Tr[A(x, x',p, ^) 7 °A(x, x',p', i/) 7 °] (A9) 
Then by simple algebraic manipulation it is trivial to show that the above trace is given by 

^L- + PzP ' z + m 2 ) r £i{AA l ) - p ± .jf ± Tv(BB')} (A10) 

Ty(AA') = 2{I v {x)I v {x')I v ,{x)I v ,{x') +I v - l {x)I v - 1 (x')I v ,- l (x)I v ,- 1 {x')) (All) 
Tr(BB') = 2{I v _ 1 {x)I v {x l )I vl - 1 {x)I v ,{x') + I v {x)I v _ l {x')I vl {x)I vl - l {x')) (A12) 

APPENDIX B 

Now to compute two-particle forward scattering matrix element for vector boson exchange case we need 

Tr{u(x,iy,p)u{x,u,p)^) (Bl) 
for direct case. It is trivial to show that the trace is given by 

^i-Tr[(A Pff7 ff + mA + Bp^-f] = ^-[p M (TrA + TrB)] (B2) 



where 



Then 



[f(x)][ 3 ^x')\ = ^l^-[iE v E v , - p z p' z )TrATrA' + P± p' ± TrBTrB'] (B3) 
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In this case the terms T1AT1B' and TrA'TiB do not contribute. In the exchange interaction for vector coupling case, 
we need 

Tr[A(x, x',p, iy)^A(x, x',p', (B4) 
By some trivial algebra, we can very easily show that this expression is given by 

-J—Tr[AA'p a p' x (^ 7 ^^) + m 2 AA'rin+VaPxBB'{ n °r 1^)] (B5) 
Then using the well known formula for the product of four 7-matrices, we have the final form of the above expression 
^^-[(p z p' z - E„E„,)Tr(AA') + 3m 2 Tr(AA') - p±.p^ ± Tr(BB')] (B6) 
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